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1. Introduction. 

Sir Michael Atiyah has made important contributions in several areas of mathematical 
physics, and these include the study of Bogomolny-Prasad-Sommerfield (BPS) monopoles; 
see, for example, Atiyah and Hitchin (1979). Static BPS monopoles are solutions of a non- 
linear elliptic partial differential equation on some three-dimensional Riemannian manifold. 
Most work on monopoles has dealt with the case when this manifold is Euclidean space 
R 3 : the equations are then completely-integrable, and can be handled by geometrical tech- 
niques. But the monopole equations on hyperbolic space H 3 are also integrable, as was 
pointed out by Atiyah (1984a,b); and in some ways, hyperbolic monopoles are simpler than 
Euclidean ones. Hyperbolic monopoles tend to Euclidean monopoles as the curvature of 
the hyperbolic space tends to zero, although this is a delicate fact which was only recently 
established (Jarvis and Norbury 1997). 

This note describes two integrable systems which are intimately related to, and were 
motivated by, hyperbolic monopoles. The first is a discrete system (or integrable mapping) , 
which in a certain sense is dual to the hyperbolic monopole system; this is the discrete 
Nahm or Braam-Austin equation. The second comes from replacing the positive-definite 
space H 3 by a Lorentzian version, namely anti-deSitter space. The Bogomolny equation 
becomes an evolution equation on this space-time, admitting soliton solutions. 

2. Hyperbolic Monopoles and the Discrete Nahm Equations. 

Motivated by the monad construction for instantons used by Atiyah et al (1978), Nahm 
(1982) discovered a kind of duality (subsequently called reciprocity: see Corrigan and 
Goddard 1984) between the monopole equations and solutions of a nonlinear ordinary 
differential equation. This ODE (described below) is called the Nahm equation. 

For monopoles on hyperbolic space H 3 of curvature — C, a variant of the Nahm con- 
struction works (Braam and Austin 1990), at least if C~ l is a positive integer. Such 
"integral" hyperbolic monopoles correspond to certain solutions of a discrete Nahm equa- 
tion: a nonlinear difference equation defined on C~ l lattice sites. This is actually a special 
case of the ADHM construction (Atiyah et al 1978) for instantons. 

One has a picture, therefore, of a correspondence (reciprocity) between monopoles 
on the 3-space of constant negative curvature — C, and solutions of the discrete Nahm 
equation on a one-dimensional lattice with lattice spacing C (provided C is the inverse 
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of an integer). The limit C — > is the continuum limit, in which the discrete Nahm 
(difference) equation becomes the Nahm (differential) equation, ft seems likely on general 
grounds that reciprocity operates, in some sense, for non-integral hyperbolic monopoles; 
but this remains an open question, fn what follows, f shall concentrate just on the discrete 
Nahm equations, and not say anything about hyperbolic monopoles. 

Let He a fixed positive integer; and let A,, Bj, Cj and Dj denote k x k matrices, 
defined for each value of the integer j. In other words, we have four k x k matrices on a one- 
dimensional lattice indexed by j. Let s± denote the forward and backward step operators 
on this lattice: so (s + ^)j = ^j+i and (s-^)j = ^j-i- For brevity, the subscript j will 
usually be omitted in what follows; thus A stands for A,-, A + = s + A stands for A, +1 , and 
so forth. 

Consider the two linear operators 

U := Cs+ + XA 

V := D_s_ - X~ X B 

(acting on a k- vector \& defined on the lattice). Here A is a constant scalar parameter, and 
the minus signs are mainly for notational convenience. The eigenvalue equations 

UV := + XA^ = CV, 

are difference equations which propagate \1/ forwards and backwards along the lattice, fn 
order for a nontrivial solution (simultaneous eigenfunction of U and V) to exist, we need 
U and V to commute, and the parameters A, [i to satisfy an algebraic relation (which 
turns out to be the vanishing of a homogeneous polynomial in these three variables) . The 
condition [U, V] = gives the discrete Nahm equation; equation (f ) is a Lax pair for it; 
and the algebraic relation defines a spectral curve (from which one can derive conserved 
quantities, show that the solutions of discrete Nahm correspond to stepping along straight 
lines on the Jacobian of this curve, etc: see Murray and Singer 1998). 
The condition that [U, V] = should hold for all A is equivalent to 

A+ = DAD' 1 , 

B+ = C- X BC, (2) 
C+D+ = DC + [A+,B+\. 
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These are the discrete Nahm (or Braam- Austin) equations. They consist of three difference 
equations for four (matrix) functions: the under-determinacy reflects the gauge freedom in 
(2). Namely, if A is a non-singular matrix on the lattice, then the system (2) is invariant 
under the gauge transformations 

A i-> AAA' 1 
B h-> ABA' 1 

(3) 

c ^ aca; 1 

D i-> A+DA" 1 . 

A gauge choice such as D = C converts (2) into a determined system: three matrix 
difference equations for three matrices. 

For completeness, let us note the correspondence with the notation of Braam and 
Austin (1990), where the equations (2) were first derived. The relation between A, B, C, 
D and their matrices /3j, 7-,- is given by A = B*, D = C* , and 

@2j = Bj , 
72j+i = Cj . 

I shall not describe the spectral curve and its consequences here; but merely list the 
eight independent conserved quantities in the k = 2 case. They are tr A, trA 2 , trS, 
tr S 2 , tr(OD), tr(ACD), tr(BCD), and tr[CD(CD - 2AB)}. Each of these expressions is 
constant on the lattice, by virtue of (2); note that they are also gauge-invariant. 

A continuum limit of (2) may be obtained as follows. Replace the integer variable j 
by t = jh, where h is the "lattice spacing", and take the limit h — > 0. Write 

C=(2h)~ 1 I+hT 3 = D, 

(4) 

B=l(T 1 + iT 2 ) = -A*, 

where the T a are antihermitian k x k matrices, / denotes the identity matrix, and star 
denotes complex conjugate transpose. Then the h — > limit of (2) is 

+ i£ a/37 [T/3, T 7 ] = 0, (5) 

which are the Nahm equations. Similarly, one obtains the standard Lax pair for the Nahm 
equations as the continuum limit of (1), if C an d A* are chosen appropriately. 



3. Reduction to a Discrete Toda System. 

It has long been known that the Nahm equation reduces to the Toda lattice. The T a take 
values in a Lie algebra, which (in the simplest case that we are considering here) is su(k). 
To reduce to Toda, one takes T3 in a Cartan subalgebra, and T\ ± 1T2 corresponding to ± 
a set of simple roots. What follows is the discrete version of this reduction. 

We express A, £?, C and D in terms of 2k lattice functions f a = f a j, Pa = Paj, where 
a = 1, 2, . . . , k, as follows: C = D = diag(/i, ^2, • • • , fk) and 
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Then the discrete Nahm equations (2) reduce to 

(Pa)+ =Pafa/fa-l , 

fa = (fa)~ + Pa+1 ~ Pa ) 

where fo is interpreted as fk, and Pk+i as pi (in other words, the index a is periodic with 
period k). The equations (7) constitute a discrete-time Toda lattice. The first example of 
such a system was that of Hirota (1977), and many other examples have been described 
more recently. 

Notice that the quantity E = ^ a =i fa ls conserved; in other words, E + = E. Let us 
define a parameter h by E = k/h 2 , and assume that f a h — > 1 as h — > 0, for each a. Then 
the ft — > limit of (7) is the differential equation 

d 2 



dt 2 



log P a =P a+ l -2p a +p ffl _ 1 , 



(8) 



which is the Toda lattice. 

Let us look in detail at the k = 2 case. Rewrite f a and p a in terms of three functions 
u, v, w, and the constant ft, according to 

fi = h~ 1 Vl + 2hw, 
f 2 = ft"Vl -2hw, 
Pi=u-v, 
p 2 = u + V. 
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Then (7) is 

w + = w + 2hu + v + 

u+ = (u- 2hwv)/y / l - 4h 2 w 2 (9) 

v + = (v — 2hwu) I a/ 1 — Ah 2 w 2 , 

which are a discrete-time version of Euler's equations for a spinning top (with an appro- 
priate choice of moments of inertia). Indeed, in the continuum limit h — > 0, (9) becomes 
dw/dt = 2uv, du/dt = —2vw, dv/dt = —2wu. These are Euler's equations, which can be 
solved in terms of elliptic functions; and (9) can be solved likewise, as follows. 
Notice that (9) admits two independent conserved quantities, namely 

= u 2 - v 2 , 

(10) 

O = 2w 2 + u 2 + v 2 - Ahwuv. 

This enables one to express u and v in terms of w (and 0, O); and then the first equation 
in (9) leads to a difference equation for w alone. Its solution is 

k 

Wj = — sn(bh) sn(bjh + c), (11) 

where k, b and c are "constants of integration" . Here k denotes the modulus of the elliptic 
functions. The conserved quantities 6 and Q are related to k and b as follows: 

= (2h 2 )- 1 [cn(bh) -dn(bh)], 
fl = (2h 2 )~ l [l - cn(bh) dn(bh)] . 

The functions u and v are then given by 

{u + v) 2 = (O - 2w 2 + A)/(l - 2hw), 

(12) 

(u - v) 2 = (O - 2w 2 - A)/(l + 2hw), 

where Aj = h~ l sn(bh) cn(bjh + c) dn(bjh + c). To see that u and v are well-defined (and 
real) for all 6, c and < k < 1, note first that 1 ± 2hw is positive. Secondly, O — 2w 2 > 0, 
from the inequality l-cndn> k 2 sn 2 . Finally, (O - 2w 2 ) 2 - A 2 = (1 - 4h 2 w 2 )0 2 > 0. As 
a final remark, note that in the limiting case k = 1, the solution is 

w = (2h)~ 1 tanh(6/i) tanh(M), 
v — u — (2/j) -1 tanh(6/i) sech(M) + tanh(6/i) tanh(6t). 
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It seems likely that this elliptic k = 2 solution corresponds to hyperbolic 2-monopoles 
with gauge group SU(2), via the Braam- Austin construction. More generally, for k > 2 one 
may speculate that discrete- Toda solutions correspond to hyperbolic /c-monopoles with Ck 
cyclic symmetry, since this is what happens for Euclidean monopoles (Sutcliffe 1996). 

4. SOLITONS IN (2 + l)-DlMENSIONAL ANTI-DESlTTER SPACE-TlME. 

Let M be a three-dimensional Riemannian manifold, with metric g and volume element r\. 
We are interested in Yang-Mills-Higgs fields on M, with gauge group SU(2) (for simplicity). 
So we have a Higgs field $ = $(x /i ) taking values in the Lie algebra su(2); here x^ = 
x 1 , x 2 ) are local coordinates on M. A gauge potential (connection) A^(x u ) determines 
the covariant derivative -D M <£> = d$/dx^ + [A^, <&]. The gauge field (curvature) is the 
su(2)-valued 2-form F^ v = D v ]. And the Bogomolny equations for $) are 

L>$ = *F, (13) 

or, in index notation, 

= \g^r} vaP F aP . (13') 

These are coupled nonlinear partial differential equations which, in general, are not com- 
pletely integrable. But (13) is an integrable system (in the sense that a Lax pair exists) 
if the metric g has constant curvature. For example, if (M, g) is Euclidean space R 3 or 
hyperbolic space H 3 , then (13) is the equation for Euclidean or hyperbolic BPS monopoles, 
respectively. 

Another possibility is for g to have Lorentzian signature — h +, and then (13) are 
evolution equations in the space-time (M, g) . Soliton solutions in the case of flat space- 
time have been studied in some detail: see Ward (1988, 1990, 1998). The aim here is to 
describe an example in curved space-time. 

There are two curved space-times with constant curvature: deSitter space with positive 
scalar curvature R, and anti-deSitter space with R < (Hawking and Ellis 1973). I shall 
deal here with the latter case only, namely anti-deSitter space (AdS). By definition, (2+1)- 
dimensional anti-deSitter space is the universal covering space of the hyperboloid H with 
equation 

U 2 + V 2 -X 2 -Y 2 = 1, (14) 
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and with metric induced from 

ds 2 = -dU 2 -dV 2 + dX 2 + dY 2 . (15) 

If, for example, we parametrize the hyperboloid Ti by 

U = seep cos 9 

V = sec p sin 9 
X = tan p cos ip 

Y = tan p sin <p 

with < p < 7r/2, then we get the metric 

ds 2 = sec 2 p(-d6 2 + dp 2 + sin 2 p dtp 2 ). (17) 

At this stage, the space-time contains closed timelike curves, because of the periodicity of 
9. Anti-deSitter space is the universal cover of Ti, in which 9 is unwound (so that 9 e R). 
Consequently, AdS, as a manifold, is the product of an open spatial disc (on which p and 
<p are polar coordinates) with time 9 G R. It is a space of constant curvature, with scalar 
curvature equal to —6. Null/spacelike infinity J consists of the timelike cylinder p = 7r/2; 
this surface is never reached by timelike geodesies. 

In what follows below, we shall also use Poincare coordinates t, x and r > 0. They 
are defined by 

t=-V/(U + X) 

r = l/(U + X) (18) 
x = Y/(U + X), 

in terms of which the metric is 

ds 2 = r~ 2 (-dt 2 + dr 2 + dx 2 ). (19) 

But these only cover a small part of AdS, corresponding to half U + X > of the hyper- 
boloid Ti. The surface r = is part of infinity X. 

The minitwistor space corresponding to AdS, or rather to the Poincare space (19), is 
CP 1 x CP 1 , which we visualize as a quadric Q in CP 3 (cf. Hitchin 1982). The points of 
space-time correspond to certain plane sections (conies) of Q. The space of all planes is a 
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CP 3 . But the relevant conies have to be real (which in this case means that their defining 
planes have real coefficients), and nondegenerate. So the space of these acceptable conies 
is the "top half" of RP S , parametrized by the real homogeneous coordinates (U, V, X, Y) 
with U 2 + V 2 - X 2 - Y 2 > 0. This RP 3 is double- covered by the original hyperboloid H, 
and is essentially the Poincare space (19). More accurately, the coordinates (t,r,x) cover 
all of RP+ except for a set of measure zero. If u> and ( are standard coordinates on the 
two CP 1 factors of Q, then the conies are 

u _ _ < ~ (™ + r 2 ) _ (Y + V)C + (X-U) 

UJ - u{c >- — — - (u + x)c + (v-Yy (20) 

where u = x + t and v = x — t. Eqn (20) expresses the correspondence between space-time 
and twistor space Q. 

The idea now is that holomorphic vector bundles V over Q (saisfying some mild 
conditions) determine multi-soliton solutions of (13) in anti-deSitter space, via the usual 
Penrose transform. In our case, the relevant vector bundles are stable bundles of rank 2, 
with Chern numbers c\ = and C2 = 2n, n being a positive integer. In the simplest case 
n = 1, the moduli space of such bundles is 5-complex-dimensional (cf. Hurtubise 1986, 
Buchdahl 1987). When we impose reality conditions, which amounts to taking the gauge 
group to be SU(2) rather than SL(2,C), the moduli space becomes 5-real-dimensional. So 
we expect, in this simplest case, to get a five-parameter family of soliton solutions, exactly 
as for the flat-space-time system (Ward 1988, 1990, 1998). This is exactly what happens 
(Hickin 1998). 

One explicit way of seeing how solutions arise is as follows: it involves a Lax pair for 
the integrable system (13). Define two operators Vi and V2 by 

Vi = rd r - 2(C - u)d u 

(21) 

V 2 = 2d v + r- 1 ((-u)d r . 

Notice that Vi and V2 both annihilate the expression (20). This is related to the fact 
that twistor space Q is the quotient of the distribution {Vi, V2} (on the four-dimensional 
correspondence space whose local coordinates are (£, r, The Lax pair involves the 

gauge-covariant version of (21), and consists of the pair of equations 

[rD r + $-2((-u)D u U = 

(22) 

[2D V + r-\C ~ u)(D r - r" 1 ^] V = 0, 
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where ip = ip(t, r,x,() is a 2 x 2 matrix. The consistency condition for this overdetermined 
system is exactly (13). 

The functions ip corresponding to n = 1 bundles can be taken to have the rational 

form 

(^-Co)p'« p (23) 

(C-Co) p-p 

where I denotes the identity 2x2 matrix, ( is a complex constant, p(t, r, x) is a row 
2-vector of linear functions of ujq = cu(Co), and p* denotes its complex conjugate transpose. 
So p has the form 

p = (auj + 6, cuj + d), (24) 

where a, 6, c, d are complex constants with ad— be ^ 0. The Yang-Mills-Higgs fields ($, A^) 
can then be read off from (22-24), and they will automatically satisfy (13). The parameters 
Co, a, b, c, d are not all significant: it is clear from (23) that an overall complex scaling of 
p will not change ip, and furthermore that multiplying p on the right by a constant SU(2) 
matrix will induce a gauge transformation on ($, A^). Removing this freedom leaves us 
with a five-real-parameter family of solutions. 

Each of these solutions represents a single soliton (lump), and the five parameters 
describe the location (2), velocity (2) and size (1) of this soliton. It is straightforward 
to write down the fields ($, A^) explicitly as rational functions of £, r, x (or U, V, X, 
Y) — but these expressions are not immediately transparent, and we shall make do with 
the following remarks. The solitons are spatially localized, in the sense that $ — > and 
F^v — > as r — > 0, i.e.. as one approaches null/spacelike infinity J. To see a simple 
example, one may set p = (a>o, 1). Then the positive-definite gauge-invariant quantity 
— tr $ 2 (which is a good one for visualizing the field) is given by 

- tr $ 2 = 8r 4 / [(r 2 + x 2 - t 2 ) 2 + 2x 2 + 2t 2 f . (25) 

The graph of this function is a single lump, with its maximum along the timelike geodesic 
x = 0, r 2 = t 2 + 1: a soliton in free fall. 

5. Concluding Remarks. 

Many of the ramifications of Atiyah's (1984b) work on hyperbolic monopoles are only now 
being addressed. In the positive-definite case, a study of the relation between hyperbolic 
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monopoles (and their symmetries) and the corresponding spectral curves is currently un- 
derway (Murray and Singer 1998). In the Lorentzian case, soliton solutions and their cor- 
responding vector bundles are being investigated (Hickin 1998); specific questions include 
multi-soliton (n > 1) dynamics, and what happens in deSitter (rather than anti-deSitter) 
space. 
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